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The  Wyngaard  and  LeMone  (1980)  model  of  interfacial 
turbulence  structure  functions  (temperature,  CT2,  and  water 
vapor,  Cq2)  in  an  entraining  mixed-layer  is  analyzed.  The 
model  indicates  that  in  the  interfacial  region  (Z  -  Z^)  Cx2 
is  proportional  to  (AX)2  Z^-2/2  9V*/AQV  where 

X  a  T  or  Q,  AX  is  the  jump  in  X  across  the  interface,  Z^  is  the 
height  of  the  interface,  and  9V*  is  the  convective  mixed-layer 
scaling  parameter  for  temperature.  Although  based  on  a  number  of 
assumptions  (referred  to  as  the  "quasi-steady"  approximation),  the 
model  is  found  to  have  more  general  application.  A  theoretical 
analysis  indicated  that  the  model  might,  not  apply  where  A0V  is 
large  (on  the  order  of  10  K) ,  particularly  for  CT2 .  A 
comparison  against  23  aircraft  profile  measurements  revealed  that 
the  model  agreed  within  a  factor  of  three. 
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I  INTRODUCTION 


>  This  report,  is  a  theoretical  and  experimental  analysis  of  a 
model  (ttyngaard  and  be M one ,  l-980^-»  used  t.o  calculate  the  refractive 
index  structure  function  parameter,  Cn2,  at.  the  inter  facial 
region  at  the  top  of  an  entraining, . turbulent  mixed  layer. 

Cn2  is  related  to  the  micrometeoro logical  structure  functions 
for  temperature,  CT^  humidity,  •  and  T-Q  covariance. 

The  mixed  layer  int.erfacial  region  is  important  for  EM 
propagation  because  iCn2  is  greatly  increased  by  large  T  and 
fluctuations  due  to  the  entrainment,  of  warm,  dry  air  from  the 
nonturbulent.  atmosphere  above  the  mixed  layer. 

Assuming  that  the  rate  of  entrainment,  is  in  equilibrium  with 
the  free  tropospheric  virtual  potential  temperature  (buoyancy) 
lapse  rate,  the  model  indicates  that.  Cx2  ia  proportional  to 
(AX) 2  0V*  Zi“,2/3/AQv  where  X  is  T  or  Q,  AX  the  jump 
at.  the  interface,  Z^  the  height,  of  the  ooundary  layer  and 
0V*  the  convective  temperature  mixed  layer  scaling  parameter. 

The  theoretical  basis  of  this  model  is  examined  and  four  data  sets 
from  the  NPS  aircraft  measurements  program  are  used  to  test,  the 
model. 
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II  THEORY 


A.  Background 

The  structure  function  parameters  for  temperature,  CT2, 
and  specific  humidity,  CQ2,  are  to  toe  evaluated  in  the 
inversion  region  by  averaging  between  heights  2  =  h0  and 
Z  =  h2  (see  Pig.  1).  The  complete  theory  was  developed  by 
Wyngaard  and  LeMone  (1980),  hereafter  referred  to  as  WL,  so  only  a 
summary  of  the  derivation  will  be  presented  in  this  report.  In  a 
few  instances  WL's  work  will  be  expanded  to  make  certain 
assumptions  and  manipulations  more  explicit. 

The  height  hQ  is  defined  as  the  top  of  the  mixed  layer 
where  w^v  *  0.  At  h2  both  fluxes  and  flux  divergences  are 
equal  to  zero.  The  average  structure  functions  are 


/  CT2  dZ 

(la) 

Jo 

r 2 

<Cq2>  =  Ah~l 

Jo  CQ2  dz 

(lb) 

where  Ah  =*  h2  -  hQ  and  the  0,  2  on  the  integral  denotes  ho,  h2  • 
The  average  structure  functions  are  related  to  their 
respective  dissipation  rates  by  the  Corrsin  equation 


<Ct2>  *  1. 6<£>“l/3<Xg> 


<Cq2 >  * 


(2a) 


1. 6<£>-i/3  <Xo> 


(2b) 


where  e  is  the  rate  of  dissipation  of  turbulent  kinetic  energy,  X 
and  Xq  are  the  scalar  dissipation  rates  (the  factor  1.6  implies 
X  is  the  rate  of  dissipation  of  temperature  variance  02). 

B.  Evaluation  of  <X> 

For  the  moment,  the  development  will  be  confined  to  the 
specific  humidity  (Q) .  The  dissipation  rate  is  calculated  from 
the  scalar  variance  budget  equation  (Q  denotes  mean  while  q 
denotes  fluctuating  specific  humidity;  later  in  the  paper  q  will 
denote  mixing  ratio,  Q/p). 

dv/dt  +  Wdv/dZ  +  d(wq)/dZ  +  2  pwq  d(Q/p)/dZ  *  -Xq  (3) 

where  v  =  q^,  W  is  the  mean  vertical  velocity  (subsidence)  and  p 
is  the  density  of  air.  Integrating  this  equation  from  hD  to 
h2»  as  in  Eq.  1,  yields 

<Xq>  =  -<Dq>  -  <Tq>  -  <PQ>  (4) 

where  D  is  the  first  two  terms  in  Eq.  3,  T  (transport)  the  third 
and  P  (gradient  production)  the  fourth.  Assuming  "quasi-steady" 
conditions,  WL  show  that  <Dq>  and  <Tq>  are  negligible  compared 
to  <PQ> ;  therefore 


<Xq>  =  -<Pq> 


(5) 
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At  this  point  the  generalized  inversion  structure  inode  1 


(Deardorff,  1979)  is  introduced 


Q  =  Q0  +  Aq  f(z); 


h0<Z<h2 


(6a) 


dQ/dZ  =  Aq  df/dZ 


(6b) 


where  f(Z)  describes  the  shape  of  the  Q  profile  in  the  inversion 
region  (assumed  to  be  the  same  for  Q  and  T)  with  f(hQ)  =  0  and 
f(h2)  -  1/  Q0  is  the  mixed  layer  value  and  Aq  the  jump  in  Q 
across  the  inversion.  Substituting  Eq.  6b  into  Eq.  5  and 
integrating  by  parts  one  obtains 


-<PQ>Ah  =  2  AQ 


d  ( wq )  /  dZ  f  dZ 


(7) 


The  mean  Q  continuity  equation 


-d (wq)  / dZ  =  dQ/dt.  +  W  dQ/dZ 


(8) 


is  used  in  Eq.  7  to  obtain 


-<P0>Ah  =  -2  AQ 


f  dQ/ dt. 

Jo 


fdZ  -  2Aq 


r 


W  dQ/dZ  f  dZ 


(9) 


The  time  derivative  of  Eq.  6a 


dQ/dt  =  aQ0/dt  +  f  dAQ/dt. 


(10) 


and  Eq.  6b  can  be  substituted  into  Eq.  la.  First  the 
"quasi-steady"  assumption  is  invoked,  setting  the  following 


conditions 


However,  since 


•  dAQ/dt  =  0 
dA0v/dt.  =  0 
dhQ/dt.  =  0 
dAh/ dt.  =  0 


dhQ/dt  =  WQ  +  Weo 


then  Eq.  11c  implies  Weo  ~  -WQ.  Assuming  constant 
divergence 

W  =  WQ  Z/hQ 


(11a 


dW/dZ  =  w0/h0 


w2  =  (1  +  a)  wc 


where  a  =  Ah/hQ  is  the  normalized  thickness  of  the  inter facial 
region.  Employing  these  relations  in  Eq.  9  and  doing  the  second 
integral  by  parts  gives 


-<Pq>  Ah  =  -2AQ  hYQ  dQ0/dt.  +  (  AQ)  2Weo(l  +  a-oZu) 


where  the  interfacial  functions  Yq  and  Zr,  are 


Despite  the  simplicity  of  Eq.  19,  WL  prefer  to  keep  the  wqG 
term  separate  in  their  development.  The  primary  reason  for  this 
is  to  simplify  the  analagous  development  for  9V  since 
w5"vo  *  Therefore,  WL  now  employ  the  "quasi-steady" 

entrainment  formula 


o 


0.8  W*  S-l/U+a) 


(20) 


where 

s  -  gr92  hg/(w*2T)  (21) 

with  Tq2  *  d0v/dZ  at  Z  =  h2  and  W*  is  the  convective 
scaling  velocity  (Z^  38  ho/0.8) 

W*3  =  g  w9vs  Zj_/T  (22) 

and  "s"  denotes  the  surface  value. 

Rather  than  make  an  explicit  substitution  for  Wqq  at  this 
point.,  one  could  keep  Weo  as  a  variable,  giving 

<Xq>  *  -2AQ  Yq  wq^/  ( ahQ )  +  (AQ)2  Weo  (l+a)SFQ/h0 

where 

Fq  -  C-2Yq  (1+a-aYQ)  +  (1+a-aZQ) ]/ ( a(l+a)S ) 

Using  the  WL  solutions  obtained  for  Eq.  24  (and  Yq  s  1/2) 

Fq  *  (6R)”1  (25) 

then 

<Xq>  »  -AQ  wq0  Yq/ ( ahQ  )  +  ( AQ )  ^^  (l+a)S/ ( 6Rh0 )  (26) 

Where  R  ■  g  A0V  hQ/(W*2T)  (21) 

The  final  result,  is  obtained  by  substituting  for  wqQ  in  Eq. 
26.  First,  the  9V  continuity  equations  at  hQ  and  at.  h2  are 


(23) 

(24) 


combined  with  the  hQ  to  h2  integral  form  similar  to  Eq.  16  to 
produce  the  relation 


d  AQy/dt.  +  Wqq  Ah(l+a)re2  "  weo  A0v(1  +  a“aYQ^  *  w^vo  (28) 
Since  WL  assume  dA0v/dt.  *  w0vo  =  0, 

A0v(l+a-ayQ)  =  F02  h(l+a)  (29) 

Using  Eq.  29  and  Eq.  18  in  Eq.  26,  one  obtains 

<XQ>  »  (AQ)2  Vleo(l+a)S(l+6-‘1)/(h0R)  (30) 

Note  that  the  first  term  in  Eq.  26  (which  was  proportional  to 
wqQ)  is  six  times  as  large  as  the  second  term  (WL  obtain  15/2 
for  this  ratio  because  they  use  two  separate  formulae  for  Weo 
which  differ  by  a  factor  of  4/5,  i.e.  6*5/4  *  15/2). 

The  development  for  temperature  is  parallel  until  the  9V 
equation  analagous  to  Eq.  26  is  reached.  Since  w0vo  =  0,  the 
final  result  is 


<X0>  =  (A0v)2(l+ct)Weo  S/(6Rh0) 
<XQ>  »  7(  AQ)  2(i+a)Weo  S/(6Rh0) 
C.  Structure  Functions 

The  final  step  in  this  process  is  to  specify 


(31a) 

(31b) 

that.  <  £>  is  one 


half  the  value  typically  found  at.  Zj_  under  convective  conditions 


<e >1/3  «  (0.2)1/3  W*  Zi”!/3 

Assuming  the  "quasi-steady"  entrainment  rate,  the  structure 
functions  become 


<Cq2>  *  3.9(AQ)2  0v*/(Zi2/3A0v) 

<CTv2>  a  0.5A9v  9v*/Zi2/3 

where  0V*  *  w0vs/w*.  The  virtual  temperature  structure 
function  is  related  to  the  temperature  structure  function, 
CT2,  by 


<Ct2>  =  2  Ti  <Ct^>/A0v 

where  Tj_  is  the  function  given  by  WL. 

One  point  worth  more  discussion  is  the  approximation 
Fq  a  (6R)_1  and  the  final  forms  of  Eq.  31.  Suppose  the 
results  of  Eq.  19  were  used  and  a  different  function  defined 

<Xq>  *  (AQ)2Weo(l+a)  S  GQ/h0 

<Xg  >  -  (A0v}2Weo(l+a)  S  FQ/h0 

where  Fq  remains  as  per  Eq.  19  but. 


(32) 

(33a) 

(33b) 

(34) 


(35a) 

(35b) 
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Gq  =  (140-  azQ)/(a  (14a  )s) 


(36) 


Using  Eq.  29  one  can  show 


FQ  -  GC  “  2YQ/R  (37) 

Following  the  calculations  WL  have  in  their  Appendix  A,  a , 
Fq  and  Gq  are  unique  function  of  R/S  (providing  the  mixed 
layer  gradient  is  zero) .  Gq  is  considerably  less  variable  than 
Fq.  The  following  formula  are  reasonable  approximations  for 
0 . 1 <r/S  <10 


a  *  0.96  R/S  -  0.11  (R/S)1*5  (38a) 

Gq  *  (1  +  0.064  \Zr7F)R-1  (38b) 

Fq  =  (1  +  0.28  \Ir/S)R~1/6  (38c) 

These  formulae  lead  to  slight  modifications  to  the  structure 
function  equations 

<CQ2>  =  3.3  (A  Q)  2‘3v*/(Zi2/3A0v)  (39a) 

<Ct2>  *  0.57  (A0V)  9v*Drp/Z  j_2/2  (39b) 

where  DT  *  1  +  0.22  ^R/S. 
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The  equations  for  Cq2  and  CT2  can  be  written  in 
various  general  forms  (now  dropping  the  bracket  notation) 


Cx2/( (AX)2DXEX)  =  1.14  0v*/(A0vZi2/3) 


or,  without  substituting  explicitly  for  Weo  and  e 


O  ,  .  O  .  0.53(1 +01  )  r.  W-n 

CX2/((AX)2DXEX)  =  ^ - 


where  D0  *  1,  En  *  3,  and 


Et  =  Ti/A0T 


D.  Discussion 

It  is  of  interest  to  ponder  the  significance  of  the  various 
"quasi-steady"  assumptions  (Eq.  11,  12,  13).  Suppose  we  exhume 
the  original  conservation  equation  integrals  from  Deardorff's  (1979) 
paper  (his  Eq.  18  and  21).  Assuming  only  horizontal  homogeneity 
and  constant  divergence,  the  general  equations  become 


Ah  dQ0/dt  +  AhYQ  dAQ/dt  -  Aq[  ( 1  -  YQ)  We2  +  YgWeo]  *  wqQ  (43a) 


Ah  d0vo/dt  +  AhYg  dA0v/dt  -  A0V[(1  -  YQ)  Wq2  +  YQ  WeQ]  =  0  (43b) 
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Thus  Eq  16  and  the  6V  analogue  can  be  reproduced  by  requiring 

We2  =  (1  +  a)  Weo  (44a) 

dAQ/dt  *  dA0v/dt  *  0  (44b) 

It  is  not  necessary  to  require  WQ  »  -W,  dh0/dt  *  dh2/dt.  »dAh/dt  = 
This  explains  why  WL  found  excellent  agreement  with  Aschurch  data 
where  W  =  0  and  Weo  ~  10  cm/s. 

Similarly,  the  general  forms  for  the  dissipation  integrals 

are 

-<Pq> Ah  *  -2 Ah  AQYQdQ0/dt  -  2AhZQAQdAQ/dt 

+(AQ)2[(i  -  Zq)  We2  +  ZqW^]  (45a) 

-<P0>Ah  -  2AhA0v  (Yq2  -  ZQ)  dA0v  /dt 

+  (A0v)2[_2YqL(1  -  Yq)  We2  +  Yq  W^]  (45b] 

+  (1  -  Zq)  We2  ♦  Zq  W^] 

which  reduce  to  Eq  14  if  the  conditions  of  Eq  44  are  met. 

Since  entrainment  and  surface  flux  tend  to  counteract  each 
other  in  the  Q  case,  it  seems  quite  reasonable  to  assume  that  the 
dAQ/dt.  and  dQ0/dt  terms  are  negligible  in  Eq  45a 

-<PQ>Ah  »  (AQ)2  [(1  -  Zq)  We2  +  ZQ  Weo]  (46) 


Instead  of  making  the  assumption  Eq  44a,  suppose  we  simply  assume 


J* 

I 

I 

I 


■  ■  i„  i i  ■  y ■  v*i».  ■  l ii  ■ '  i> ii  v 


we2  =  weo  "  aw6vg/A0, 


(47) 


which  is  the  standard  cloud-free  result  from  Lilly  (1968)  where 
typically  a  »  0.2.  Then  one  can  easily  show  that 


<Cq2>  -  3.3  (AQ)2  0v*/(2i2/3A0v) 


(48) 


which  is  identical  to  the  WL  result  as  expressed  in  Eq.  39a 1  In 
other  words,  the  combination  of  "quasi-steady"  assumptions 
We2  ■  (1  +  a )weo  aad  we2  "  w6vs/(r0ho)  are 

equivalent  to  the  assumptions  of  Eq.  47  even  though  they  may  imply 
vastly  different  entrainment  rates. 

If  one  uses  the  assumptions  of  Eq.  47  and  parallels  the  WL 
development,  then  the  equivalent  to  Eq.  18  is 


wcIo  *  ~  Weo  (49) 


and  the  equilibrium  condition  form  the  9V  equation  is 


vo  *  weo  ~  ^®v^  (50) 


which,  assuming  wQ vo  ■  0,  gives 


a  -  R/S  (51) 


The  results  for  9V  are  also  interesting  oecause  it  is  not 
clear  that  the  dA0v/dt  term  should  be  negligible  compared  to  the 
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other  terms  in  Eq  46  b.  Suppose  we  let 


-<Pg>  Ah  =  A  +  B 


(52) 


Then  the  dA0v/dt.  term  is  small  if  A/B  is  small  (returning  to  the 
"quasi-steacfy"  format) 


A/B 


h0(YQ2  -  Zq)  dA0v/dt  6(R/S) 
A0V  (1  +  a)  weo 


Since  Yq2  -  Zq  ~  -0.1,  we  can  write 


(53) 


A/B 


-0.6  h0(R/S)  dA0v/dt. 
(1  &)  W  go  ^®v 


(54) 


The  magnitude  of  A/B  can  be  examined  by  using  the  general 
relation 


dA0v/dt.  *  -d0o/dt  +  rQWe2 


(55) 


and  writing  a  simple  entrainment  formula  (e.g.  "quasi-steady”) 


We2  *  w 0 v s /  ( ^  g ^o  ^ 


(56) 


The  integral  of  the  conservation  equation  from  Z  *  0  to  2  *  hQ  gives 


(57) 


ho  dQvo/ “  w6vs  + 

therefore 

dA0v  /dt.  »  -  Weo  A0v/ho 

using  Eq.  54  we  find 

0.6 

A/B  =  - - -  R/S  (59) 

( 1+ct ) 

A  good  example  is  the  Aschurch  data  quoted  by  WL  where  Eq.  57 
was  shown  to  be  applicable.  Since  R/S  =  0.3  for  that  data,  A/B  * 
0.15  and  dA0v/dt  is  negligible. 

Certainly  the  conditions  set  by  WL  are  consistent  with 
neglecting  dA0v/dt.  It  is  not  clear  how  to  identify  conditions 
where  this  assumption  is  invalid.  Eq.  54  cannot  provide  much 
guidance  because  it  is  based  on  solutions  to  Eq.  28  with 
dA0v/dt  *  0.  It  is  interesting  that  in  the  conditions  where  the 
WL  equations  for  "quasi-steady'’  entrainment  are  expected  to 
breakdown  (A0V  large,  R/S  >1)  then  the  Lilly  type  relations  give 
the  same  results  for  Cq2.  if  the  dA0v/dt  term  becomes 
important,  then  one  anticipates  the  WL  formulation  will 
underestimate  Ot>2. 


Ill  ATMOSPHERIC  DATA 


A.  Measurement.  Techniques 

The  measurements  were  made  using  a  single  engine  Bellanca 
Viking  aircraft  operated  by  Airborne  Research  Associates  of 
Weston,  MA.  The  instrumentation  and  data  processing  have  been 
previously  described  in  detail  (Fairall  et.  al.,  1980;  Schacher 
et.  al. ,  1980)  so  only  a  brief  summary  is  given  here. 

i)  Mean  temperature,  T:  platinum  resistance  sensor  with 
standard  aircraft  mount. 

ii)  Mean  humidity,  Q:  cooled  mirror  dew  cell. 

iii)  Mean  windspeed,  U:  estimated  at  the  surface  from  the 
sea  state  and  DMV  navigational  aid.  The  present  LORAN  system 
was  not  available. 

iv)  Sea  surface  temperature,  Tg:  Barnes  PRT-5  IR 
radiometer. 

v)  CT2:  microthermal  sensors  (4.5  um  dia.  tungsten) 
in  the  paired  configuration. 

vi)  Cq2:  Lyman-alpha  fast  humidiometers  using  the 
inertial  subrange  filter  method.  Absolute  calibration  based 
on  comparison  with  a  microwave  refractometer . 

vii)  e:  hot  wire  (4.5  ym  dia.  tungsten)  constant 
temperature  anemometer.  The  inertial  subrange  filter  method 
was  used. 

B.  Surface  Fluxes  and  Turbulence  Scaling  Parameters 

Surface  fluxes  were  evaluated  from  aircraft  measurements 
using  two  methods;  a)  bulk  aerodynamic  and  b )  dissipation 
(inertial  subrange).  The  fluxes  are  defined  in  terms  of  the 
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following  scaling  parameters: 


momentum: 


puws  =  -pu*2 


(60a) 


sensible  heat:  pCp  w9s  =  -pCp  u*T* 


moisture: 


pqws  =  -pu*q* 


(60b) 


(60c) 


The  momentum  flux  is  also  referred  to  as  the  Reynolds  stress,  x. 
Note:  the  bulk  method  was  not  used  overland. 

1.  Bulk  aerodynamic  method. 

The  exact  details  were  described  in  a  recent  paper  (Davidson 
et  al,  1981).  Using  Eq.  4a  from  that  paper,  one  can  relate  the 
values  of  some  meteorological  variable  (temperature,  moisture  or 
wind  speed)  at  the  sea  surfce,  Xs,  and  at  some  height  Z  in  the 
surface  layer,  Xz,  to  the  scaling  parameter,  X*: 


u*  -  uzk[Jtn  (Z/Z0)  -  fu  (Z/L)]“l 


T*  *  (Tz  -  Tg)  0k[*n  (Z/ZoT)  -  f-r  (Z/L)]-l 


q*  =  (qz  -  qa)  3*  C  In  (Z/ZoT)  -  %  (Z/L)]-1 


(61a) 


(61b) 


(61c) 


where  ZQ  and  Z0>i<  are  roughness  lengths,  L  is  the  Monin-Obukhov 
length,  3  and  k  are  constants,  and  fu  and  %  are  empirical 


functions . 


2.  Dissipation  method. 

The  dissipation  method  relies 
of  inertial  subrange  turbulence  to 
parameters  (Fairall  et  a_l.  ,  1980). 


on  semi-empirical  relationships 
surface-layer  scaling 
The  equations  are 


u*  =  [(6k  Z  )/<j>  (Z/L)]1/3 


(62a) 


T*  =  [Z2/3  CT2/f (Z/L) 31/2 


(62b) 


Q*  -  [Z2/3  Cq2/ ( A  f(Z/L)]l/2 


(62c) 


where  e  is  the  dissipation  rate,  $  and  f  are  empirical  functions, 
and  A  is  a  constant.  Because  the  structure-function  parameters 
C«p2  and  Cq2  are  related  to  the  square  of  the  scaling 
parameter,  a  sign  ambiguity  exists.  This  can  usually  be 
eliminated  by  examining  the  low-level  height  dependence  of  £, 

Cq2  and  CT2  because  the  functions  <j>  and  f  have 
characteristic  profiles  for  stable  and  unstable  conditions. 

3ot.h  methods  yield  accuracies  on  the  order  of  10%  for  u*, 
+0.02°C  for  T*  and  +0.02  g/'m3  for  Q*  (note:  q*  =  Q*/p). 

C .  Data  Sets 

The  data  given  in  this  report,  were  oct.ained  in  four  field 
programs: 

i)  Panama  City  (PC),  1978  (more  detail  available  in 
Fairall,  1979)  over  the  Gulf  of  Mexico  m  Florida. 

ii)  White  Sands  (WS),  1979.  Two  profiles  over  the  desert, 
under  highly  convective  daytime  conditions. 


IS 


iii)  MAGAT  (MG),  1980  (more  detail  available  in  Fairall, 


1980)  in  the  Monterey  Bay  area. 

iv)  Bahamas  (BH),  1980.  A  series  of  profiles  taken  near 
Andros  Island  in  the  late  fall. 

The  complete  data  sets  were  examined  to  remove  profiles  that 
encountered  boundary-layer  clouds.  A  total  of  23  profiles  were 
selected.  Graphs  of  the  mean  and  turbulence  profiles  for  each 
case  are  given  in  Appendix  A.  A  summary  of  the  basic  scaling 
parameters  is  given  in  Table  1. 
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IV.  RESULTS 


A  summary  of  the  secondary  scaling  parameters  used  for  the 
NPS  data  set  is  given  in  Table  2.  Also  shown  in  Table  2  is  a 
comparison  of  the  measured  and  model  assumed  values  for  t  at  the 
inversion.  With  very  few  exceptions,  the  model  assumption  (Eq. 
32)  is  very  good.  The  entrainment  velocities  calculated  from  the 
"quasi-steady"  assumption  used  by  WL  (Eq.  20)  and  the  more 
conventional  parameterization  of  Lilly  (1968). 

Weo/W*  =0.2  ev*/A0v  (63) 


are  also  calculated. 

In  Table  3  are  the  measured  values  of  CT2  and  Cq2  at 
the  inversion  plus  their  normalized  forms 

XX  *  Zi2/3  CX2/((AX)2  Dx  FX)  (64) 

taken  from  Eq.  40.  According  to  WL  (Eq.  26),  the  theoretical 
value  is 


Ic  *  1.14  9V*/A0V  (65) 

which  is  the  same  for  T  and  Q. 

A  direct  comparison  of  measured  and  calculated  values  of 
Crp2  and  Cq2  is  given  in  Fig.  2.  The  model  predicts  the 
measurements  within  a  factor  of  three.  The  uncertainty  is 
slightly  greater  than  the  factor  of  two  suggested  by  WL  but 
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Measured  values  of  the  interfacial  structure  functions 
(CT2  and  CQ2)  and  their  resultant  values  for 

"X  “  Zj.2/2  CX2/((AX)2  Dx  F x^ where  X^T  . 

or  u  'l'nese  are  compared  with  theoretical  values,  Ic# 

using  the  "steady"  and  "Lilly"  entrainment  values. 
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Schematic  representation  of  the  convective  boundary 
layer  (taken  from  Wynyaara  an<i  LeMoiie,  1980)  with  . 
inrerfacial  Layer  snow i  rm  h  ,  h- ,  ih ,  l  Loves 

and  jumps.  ‘ote  that  l-.\y  -  ^  v  ( h  2  )  -  .''v(h,J) 
is  positive  while  i>_i  is  negative. 


includes  various  measurement  errors  and  uncertainties.  Note  that 
the  Cq2  data  has  a  greater  range  of  values  tnan  .  This 
is  consistent  with  the  WL  model.  If  we  examine  the  function 

H  -  Zi2/3  Cx2/(dx  9v*)  (66) 

then 

Ht  =  FT  A0  v  (67a 

Hq  =  (AQ)  2/A0  v  (67b 

A  graph  of  HT  and  Hq  is  shown  in  Fig.  3  for  a  typical  range  of 
A0V  and  AQ  from  the  NPS  data  set.  Note  that  HT  varies  roughly 
from  2  to  9  while  Hq  varies  from  4  to  72. 

The  entrainment  parameterization  was  tested  (Fig.  4)  by 
plotting  measured  values  of  Ix(Eq.  65)  against  the  model  value 
(Eq.  66)  which  is  based  on  the  entrainment  formula  given  by  WL 
(Eq.  26).  This  plot  gives  a  much  higher  correlation  than  a 
similar  graph  (not  shown)  using  the  more  traditional  formula  due 
to  Lilly  (1968),  Eq.  62,  which  gives 

Ic'  (Lilly)  =  0.18  (l+o)r0  ZL  9v*/(A0v)2  (68) 

This  is  not  really  significant  because,  when  used  in  proper  com¬ 
bination  with  Eq.  48,  the  Lilly  formulation  also  leads  to  Eq.  66. 

In  order  to  look  for  systematic  errors,  the  ratios  (Rt  and 
Rq)  of  measured  to  model  values  of  C ^  and  Cg^  were 
calculated  and  plotted  against  .10  v  (Fig.  5).  A  simple 
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Theorer  ioa  L  uxoresaiuii  cor  H^.  ami  H,  (dq.  of.) 
illustrating  r  ne  difterenoe  between  the  dependence  ;>£ 
Cr^.2  ana  on  A0V  ana 
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Iy .THEORY 


IQ, THEORY 


Rt  (Meas./ Theory! 


The  measured  value  or  divided 

value  as  a  function  of 


Rq  (Meas./Theoryl 


log-average  yields  RT  =  1.15  and  Rq  =  1.3.  Figure  5a  weakly 
suggests  that  the  model  underestimates  C<j>2  (large  ftp)  when 
A0V  is  small  while  it  overestimates  when  A0V  is  large  (the 
Cq2  data  is  too  sparse  to  clear  up  this  question) .  This  could 
be  due  to  an  error  in  the  estimation  of  A0V  and  AQ  (admittedly 
rather  subjective).  An  examination  of  Fig.  3  suggests  that  a 
reasonable  adjustment  of  A0V  (several  tenths  K)  will  not  move 
the  data  points  substantially  closer  to  the  Rrp  =  1  midline. 
Another  possibility  is  that  Eq.  20  tends  to  overestimate  Weo 
when  A0V  is  large  while  underestimating  for  small  A0V. 

Given  the  considerable  scatter  in  the  results,  the 
uncertainties  in  the  estimation  of  A0V  and  AQ  from  measured 
profiles  and  the  insensitivity  of  Cp2  to  A0V  and  AQ  it  is 
suggested  that  a  simplified  formula  for  CT2  can  be  used  for 
application  to  radiosonde  quality  data.  If  one  assumes  (based  in 
Fig.  3)  that  HT  -  5,  then 

cT2  2  5  ev*  Zi“2/3  (69) 

Based  on  the  NPS  data  set  this  approximation  appears  to  be  at 
least  as  accurate  as  the  more  complicated  formula  (Fig.  6). 


V  CONCLUSIONS 


The  Wyngaard-LeMone  inversion  layer  scaling  has  been  examined 
theoretically  and  tested  against  a  data  set  obtained  by  NPS 
investigators  in  cooperation  with  Airborne  Research  Associates. 

The  theoretical  examination  indicated  the  following: 

i)  The  WL  theory  is  more  general  than  is  implied  by  the 
strict  assumptions  of  the  "quasi-steady"  theory. 

ii)  The  WL  development  can  be  simplified  slightly,  leading 
to  modest  adjustments  of  the  normalization  constants. 

iii)  The  steady  state  assumption  that  dAQ/dt.  is  negligible  is 
reasonable  under  most  conditions.  The  assumption  that  dA0v/dt. 

is  negligible  may  not  be  justified  when  R/S  >  1. 

The  examination  of  the  atmospheric  data  indicated  the 
following: 

i)  The  assumption  that  e  at  the  inversion  is  proportional 
to  a  fixed  fraction  of  the  surface  buoyancy  flux  was  quite 
reasonable . 

ii)  The  WL  model  predicted  the  measured  value  of  C ^  and 
Cq2  to  within  a  factor  of  three. 

iii)  Some  evidence,  though  statistically  weak,  was  found  to 
suggest  the  model  overestimates  the  structure  functions  for  large 
A0V  ( >  3K) )  while  it  underestimates  for  small  A0V  ( <  2K) .  On 
the  other  hand,  this  could  be  a  manifestation  of  the  Stein  effect 
for  comparison  of  data  sets  subject  to  error  where  small  values 
are  usually  overestimated  and  large  quantities  are  usually 


underestimated. 


Based  on  these  results,  it  is  obvious  that  a  major  weakness 
of  the  model  is  its  reliance  on  an  entrainment  formulation  that  is 
too  restrictive.  The  two  extremes  of  the  buoyancy  jump  (A0V) 
may  involve  different  entrainment  regimes  (e.g.  encroachment, 
convective  instability  or  the  Lilly  formulation).  It  would  also 
be  useful  to  include  the  effect  of  inversion  windshear  on  We  and 
on  the  structure  functions.  Another  area  of  investigation  might 
be  stable  surface  layers.  These  may  be  very  important  for  surface 
optical  propagation  because  C-p2  values  are  often  sizeaDle  and 
is  usually  small  (on  the  order  of  100m) . 
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APPENDIX  A 


This  appendix  contains  graphs  of  mean  (Qv,  q)  and 
turbulence  (C«p2 ,  Cq^,  e)  profiles  for  each  of  23  data 
sets.  The  site  designations  are  defined  in  Section  III-C.  The 
abstraction  of  this  data  to  obtain  the  relevant  parameters  (Tables 
1,  2,  3  in  the  main  text)  is  described  in  Section  III. 
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Figure  A6b.  Turbulence  profile  lor  PC  12/10  1523. 
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